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Correlations of non-photonic electrons with charged hadrons are studied using the 
PYTHIA Monte Carlo event generator in p+p collisions at y^sjvjv = 200 GeV. We 
propose experimental methods to estimate the relative contributions to non-photonic 
electrons from charm and bottom decays. We present the measurement of azimuthal 
correlations between non-photonic electrons and charged hadrons in p+p collisions at 
ytspffj = 200 GeV from STAR. The results are compared to PYTHIA simulations to es- 
timate the relative contributions of D and B meson semi-leptonic decays to non-photonic 
electrons. Our preliminary result indicates that at ~ 4.0 — 6.0 GeV/c the measured B 
contribution to non-photonic electrons is comparable to D contribution, and is consistent 
with the FONLL calculation. 



1. Introduction 

High transverse momentum particles in central Au+Au collisions are found to be 
significantly suppressed compared to that from p+p collisions This is consistent 
with the picture that high energy partons are quenched when propagating through 
the dense medium created in central nucleus-nucleus collisions^ and this suppres- 
sion is thus a probe for the Quark Gluon Plasma. Light hadron sup pres sion at 
RHIC has been described by a gluon radiation energy loss mechanismF*^ Recent 
RHIC datcP^ show that the suppression of non-photonic electrons (electrons from 
heavy quark decays) has almost the same magnitude as that of light hadrons in 
central Au+Au collisions. This implies that heavy quarks may lose a substantial 
amount of energy, which cannot be explained by current theoretical predictions 

|Q|Q| 

based on gluon radiation as the dominant mechanism for energy loss. 121 ^ 1 Heavy 
quark would lose much less energy because the gluon radiation is suppressed at 
small angles due to the mass of heavy quarks (dead cone efFectP^J). Recent calcula- 
tions including b oth glu on radiation and collisional energy loss increase the heavy 
quark energy loss ff^^ In both scenarios electrons from charm decays are always 
more suppressed than those from bottom. However, theoretical calculations contain 
large uncertainties with respect to the relative contributions to non-photonic elec- 
trons from charm and bottom quarks. Another important tool to study the Quark 
Gluon Plasma properties is the angular anisotropy V2 ■ Observation of heavy quark 
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hydrodynamic flow is a unique probe for the early stage of a partonic phaseff^ A 
nonzero non-photonic electron v 2 has been measured for p? < 2.0 GeV/c, while 
at higher Pr (although uncertainties are still large) the v-i is observed to decrease 
with pr^^ The same D and B hadron v% can lead to very different non-photonic 
electron v-i due to the different decay kinematics between D and B hadronsP^ 
The quantitative understanding of the non-photonic electron suppression at high 
Pt and V2 requires the knowledge of the relative charm and bottom contributions 
to non-photonic electrons. 

In this paper, correlations of non- photonic electrons with charged hadrons are 
studied using the PYTHIA event generator in p+p collisions at ^/snn — 200 GeV. 
We develop an innovative method which uses the azimuthal correlations between 
non-photonic electrons and charged hadrons to estimate the relative D and B con- 
tributions to non-photonic electrons. We present the preliminary results on the 
measurement of azimuthal correlations between high-p^ non-photonic electrons and 
charged hadrons in p+p collisions at 200 GeV from STAR. We will use comparisons 
of the experimental results with PYTHIA simulations to estimate B and D decay 
contributions to non-photonic electrons as a function of pr for pt > 2.5 GeV/c. 



2. PYTHIA Simulation 

We have studied the correlations between non-photonic electrons and inclusive 
charged hadrons in p+p collsions at ^/s^n = 200 GeV from PYTHIA i;6.22pH 
The PYTHIA parameters have been tuned in order to simultaneously describe the 
STAR measurements of the D meson pt and of the single non-photonic 

electron pt distributional The parameters for charm quarks are: PARP(67) = 4 
(factor multiplied to Q 2 ) , < k t > = 1.5 GeV/c, m c = 1.25 GeV/c 2 , Kf actor = 
3.5, MSTP(33) = 1 (inclusion of K factors), MSTP(32) = 4 (Q 2 scale), CTEQ5L 
PDF and S fragmentation function. The parameters for bottom quarks are the same 
as for charm quarks except mb = 4.8 GeV/c 2 . The details of PYTHIA parameter 
setup are discussed in Reference!^ 

The left panel of Fig. [1] shows the A(p distributions between non-photonic elec- 
trons and inclusive charged hadrons with six electron trigger pt cuts and the as- 
sociated hadron pr > 0.3 GeV/c. The distributions are scaled by the number of 
electron triggers. The solid curves and dashed curves are for B decays and for D 
decays, respectively. There is a significant difference between B and D meson de- 
cays in the near-side correlations. The width of near-side peak for electrons from 
D decays is much narrower than those for the B decays. The wide width from B 
meson decays is due to the larger energy release (Q value) in the B meson semi- 
leptonic decays leading to a broad angular correlation between daughter hadrons 
and electrons in the laboratory frame. For an electron at high px from B meson 
decays, the B meson does not have to be at high transverse momentum because 
the electron can get large momentum from the 6-quark mass. In the case of D 
meson decays, the D meson needs to have a large momentum in order to boost 
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Fig. 1. (color online) Left: PYTHIA simulations of angular correlations between non-photonic 
electrons and charged hadrons with six electron trigger px cuts and associated hadron pj- > 0.3 
GeV/c. Solid lines show electrons from B meson decays. Dashed lines show electron from D 
meson decays. Right: PYTHIA simulations of summed px distributions of charged hadrons around 
triggered non-photonic electrons with three trigger pt cuts. Solid lines show electrons from B 
meson decays. Dashed lines show electrons from D decays. 




Fig. 2. (color online) The distributions of electron pt versus its parent pt (left for D decays and 
right for B decays). 

the daughter electron to a high px- Fig. [5] shows the distributions of electron px 
versus its parent pr- We found the difference in the near-side correlations between 
D decays and B decays is largely due to the decay kinematics, not the production 
dynamics. Variations on fragmentation function from default Peterson function to 
8 function and on the parameter PARP(67) from 1 to 4 do not change the shape 
of the correlations in a significant way. This difference can help us estimate the 
relative B and D contributions to the yields for non-photonic electrons. We will 
use the A<p distributions from B decays and D decays in PYTHIA to fit the Aip 
distributions from the real data and thus measure the B contribution, B j (B + D) . 

We further studied the particle production within a cone around triggered high 
Pt electrons from heavy quark decays. We focused on the scalar summed px dis- 
tributions of inclusive charged hadrons in the cone (pt refers to the transverse 
momentum in the laboratory frame). Here the cone is defined by \rjh — Ve\ < 0.35 
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Fig. 3. (color online) (a)Invariant mass distributions of e + e — pairs (OppSign, grey filled area) and 
combinatorial background (SameSign, red solid line) in p+p collisions. (b)The ratio of inclusive 
electron to photonic background as a function of in p+p collisions. 



2.5 < P T (trig) < 3.5 GeV/c, P (asso) > 0.3 GeV/c 



3.5 < P T (trig) < 4.5 GeV/c, P T (asso) > 0.3 GeV/c 



0.1 

0.05 




: STAR Preliminary 


X 2 / ndf 
B/(B+D) 

i 


29.09 / 29 
0.2819 + 0.0649 


Ik 


r^jc — . 

1^ 




- 1 .... i .... i . 


i . . . . 





Z 2 'ndf 41.45/29 
B/(B+D) 0.4087 + 0.0887 




a? 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 




X 2 'ndf 51.2/29 
STARPrchmmar, B;(B+D) . 3865 + 00630 






A9 



5.5 < P T (trig) < 6.5 GeV/c, P T (asso) > 0.3 GeV/c 
0.45i — 



20.11 /29 
0.4654 + 0.0995 




Fig. 4. (color online) Aifi non _ p f lo distributions and the comparison to PYTHIA simulations for 
four electron trigger cuts with associated hadron pT(assoc) > 0.3 GeV/c. The data are shown as 
dots, the simulations are depicted by blue dashed lines for B meson decays and red dotted lines 
for D meson decays. The green solid curves are the fits to data points using PYTHIA curves. 



and \iph — ¥> e | < 0.35 (77 is pseudo-rapidity and ip is azimuthal angle). The summed 
Pt distributions of inclusive charged hadrons in three triggered electron px ranges 
are shown in the right panel of Fig. [TJ The distributions are scaled to unity. The 
dashed lines are for D decays and the solid lines are for B decays. We also can see 
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that there is a significant difference between B decays and D decays. The summed 
Pt distributions for D meson decays are much wider than those for B meson decays. 
This difference can also be used to distinguish B and D decay contributions. It will 
be more valuable for the small acceptance experiments to investigate the B decay 
contribution using summed pt distributions. 

3. STAR Data Analysis 

The data used in this analysis is p+p events at ^snn — 200 GeV recorded by 
STAR in RUN V. The main detectors utilized in this analysis are the STAR Time 
Projection Chamber (TPC)^- and the STAR Barrel Electromagnetic Calorime- 
ter (BEMC) with the Shower Maximum Detector (SMDJ2H. To obtain sufficient 
statistics at high-pT, we used high tower triggers corresponding to an energy de- 
position of at least 2.6 GeV (HT1) and 3.5 GeV (HT2) in a single tower of the 
BEMC. Around 2.4 million HT1 events and 1.7 million HT2 events were used for 

Electron identification was car- 
ried out by combining ionization 
energy loss in the TPC with energy 
deposition in the EMC and shower 
profile in the SMD. The measure- 
ment of the ionization energy loss, 
dE/dx, for charged tracks in the 
TPC gas is used to identify elec- 
trons in the first stage. Requiring 
the dE/dx values of the selected 
tracks to be near the expected elec- 
tron band in the region px > 2.0 
GeV/c rejects a significant fraction 
„. K , , !• wu d/jdi n\ »' f of the hadron background. After 

b lg. 5. (color online) 1 he Bj (B+JJ) ratio as a mnction ° 

of the fit range in Atp. Scale factors are used in order extrapolating the TPC tracks to 

to separate the data points well. the BEMC, we require the ratio 

p/E to be less than 1.5 using the 
momentum information from the TPC, and the tower energy information from 
the BEMC, E. Electrons will deposit almost all of their energy in the BEMC while 
this is not true for hadrons. Further hadron rejection is provided by the shower max- 
imum detector (SMD), which allows us to cut on the shower size with high spatial 
resolution. We require the profile of the electro-magnetic shower to be within the 
expectation for electrons. Combining the power of TPC and BEMC, we can achieve 
an inclusive electron sample with purity > 98% in the pt region up to 6.5 GeV/c. 
More detailed descriptions of the electron reconstruction techniques can be found 
in References. ^221 

The dominant photonic electron background is from photon conversions and 
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Fig. 6. The fit \ 2 as a function of the B/(B + D) ratio. 



7T°, rj Dalitz decays, whose e + e~ pairs have small invariant masses. We combine 
the electron candidates with tracks passing a very loose cut on dE/dx around 
the electron band. The invariant mass distribution ol e + e~~ pairs (OppSign) is 
depicted by the grey filled area in panel (a) of Fig. [3J The OppSign contains the 
true photonic background as well as the combinatorial background, where non- 
photonic electrons may be falsely identified as photonic electrons. The combinatorial 
background, which is small in p+p collisions, can be estimated by calculating the 
invariant mass of same-sign electron pairs (SameSign) shown as red solid line in 
panel (a) of Fig. [3J A cut of mass < 0.1 GeV/c 2 rejects most of photonic background. 
Panel (b) of Fig. [3J shows the ratio of inclusive electron to photonic background 
as a function of px- The bars (boxes) on the data indicate the size of statistical 
(systematic) errors. A significant excess of electrons with respect to the background 
has been observed, where the excess electrons are mostly from heavy quark semi- 
leptonic decays. In RUN V there was a larger amount of material producing photon 
conversions in the STAR experimental configuration. This led to inclusive electron 
to photonic background ratios from RUN V systematically lower than those from 
RUN III (IV) P 

In order to extract the angular correlation between non-photonic electrons and 
charged hadrons, we start with the semi-inclusive electron sample. We remove the 
OppSign sample after the mass cut from the inclusive electron sample. The re- 
maining electrons form the " semi- inclusive" electron sample. The relationship of 
these samples is: semi-inc = inc - OppSign with the mass cut — inc — (reco-pho 
+ combinatorics) — inc — (pho — not-reco-pho + combinatorics) = non-pho + 
not-reco-pho — combinatorics. Therefore the signal can be obtained by the equa- 
tion. ^ip non ^pf lo — A(p scrai~inc ^^Pnot — reco~pho~ > x~ ^^combinatorics- ^^Pnot — reco—pho 



can be calculated using A(p reco - p h by an efficiency correction after removing the 
photonic partner of the reconstructed-photonic electron. The reason to remove the 
photonic partner is that for the reconstructed-photonic electron the photonic part- 
ner is found while for not-reconstructed-photonic electron the partner is missing. 
The resulting e-h correlations for reconstructed photonic electrons and not recon- 
structed photonic electrons cannot be related to each other by an efficiency cor- 
rection factor alone. Therefore Atp not - re co-pho can be obtained by the equation: 
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A^not-reco-pho = (1/e - 1) ^reco-pho-no-partner, where e is the photonic electron 
reconstruction efficiency and Aip reco - p h -no-partner means reconstructed photonic 
electron azimuthal correlations with charged hadrons after removing the photonic 
partners. The corresponding efficiency can be calculated from simulations and is 
~ 70%. 

4. Results and Discussions 

Fig. [4] shows the Atp non _ pho dis- 
tributions and the comparison to 
PYTHIA simulations for four elec- 
tron trigger cuts with associated 
hadron pt(cissoc) > 0.3 GeV/c. 
The data are shown as dots, the 
blue dashed curves and the red dot- 
ted curves are from PYTHIA sim- 
ulations for B decays and for D de- 
cays, respectively. We use PYTHIA 
v 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 curves to fit the data points with 
P T (GeV/C) the B contribution as a parame- 

„. _ , . .. . _ ., ter in the fit function: Aip exv — 

big. 7. (color online) the relative B contribution to 

non-photonic electrons as a function of electron pt- R X A(fg + (1 R) X Aipjj, where R 

is the B/{B + D) ratio. The fits are 
shown as green solid curves in Fig. [4] The fit results are consistent within statistical 
errors when we vary the fit range in Aip from ±1 to ±7r as shown in Fig. [5l As 
a cross check, we fixed the B/(B + D) ratio to see how the fit x 2 changes. Fig. [6] 
shows the fit % 2 as a function of the B / (B + D) ratio with the fit range in Aip from 
— 1.5 to 1.5 for two trigger pr cuts. The x 2 1S sensitive to the Bj(B + D) ratio, and 
so is it for other trigger px cuts and for other fit ranges. As a systematic check, 
we allowed an overall normalization factor in the fit function to float. It gives a 
normalization factor close to unity, and consistent B / (B + D) ratios. 

The B semi-leptonic decay contribution to non-photonic electrons as a function 
of pt is shown in Fig. [7] The bars show the size of statistical errors. The grey 
band indicates the data uncertainties including statistical errors, and systematic 
uncertainties introduced by photonic electron reconstruction efficiency (dominant), 
different fit ranges and different fit functions. The blue solid curves show the range of 
relative bottom contribution from the Fixed Order Next-to-Leading Log (FONLL) 
calculations.^ The dashed line is the B/(B + D) ratio in FONLL used for the 
default non-photonic electron Raa calculation. The preliminary STAR data are 
consistent with the FONLL calculation. Taking the observed non-photonic elec- 
tron suppression in central Au+Au collisions at RHIC together with the measured 
B/(B + D) ratio in p+p collisions suggests that B mesons may be suppressed 
in dense medium. A calculation from Adil and Vitev predicts suppression of B 
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mesons comparable to that of D mesons at transverse momentum as low as px ~ 10 
GeV/cP 

5. Conclusion 

In conclusion, correlations of non-photonic electrons with charged hadrons are stud- 
ied using the PYTHIA model in 200 GeV p+p collisions. The measurement of az- 
imuthal correlations between non-photonic electrons and charged hadrons in p+p 
collisions at ^/snn — 200 GeV from STAR has been presented. The experimental 
results are compared with PYTHIA simulations to measure the relative B contribu- 
tion to non-photonic electrons. Within the present statistical and systematic errors, 
the preliminary data analysis based on PYTHIA model indicates at pr ~ 4.0 — 6.0 
GeV/c the measured B contribution to non- photonic electrons is comparable to 
the D contribution and that the measured B/(B + D) ratios are consistent with 
the FONLL theoretical calculation. Together with the observed suppression of non- 
photonic electrons in central Au+Au collisions, the measured Bj(B + D) ratios 
imply that the bottom quark may be suppressed in central Au+Au collisions at 
RHIC. 

This work is partially supported by NSFC under the Grant No. 10575042, 10610285. 
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